Soft and collinear parton radiation in heavy quark production^ 
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When the energy of the heavy quark is comparable with its mass, it is natural to attribute this 
heavy quark to the hard part of the reaction. At large energies, this approach is impractical due to 
£SJ ' large logarithms from intensive QCD radiation affecting both inclusive and differential observables. 

We present a formalism for all-order summation of such logarithms and reliable description of heavy- 
quark distributions at all energies. As an illustration, we calculate angular distributions of B mesons 
produced in neutral-current events at large momentum transfers at the ep collider HERA. 

In recent years, significant attention was dedicated to exploring properties of heavy-flavor hadrons pro- 
duced in lepton-nucleon deep inelastic scattering (DIS). On the experimental side, the Hadron-Electron Ring 
Accelerator at DESY (HERA) has generated a large amount of data on the production of charm and bottom 
mesons in neutral current interactions (cf. jjj and references therein). At present energies (about 300 GeV 
in the ep cm. frame), a substantial charm production cross section is observed in a wide range of Bjorken 
x and photon virtualities Q 2 , and charm quarks contribute up to 30% to the DIS structure functions. 

On the theory side, Perturbative Quantum Chromodynamics (PQCD) provides a natural framework for 
the description of the heavy-flavor production. Due to large masses M of charm and bottom quarks (M 2 ^> 
t"~- Aq CD ), the renormalization scale can be always chosen in a region where the effective strong coupling as is 

small. Despite the smallness of as , perturbative calculations in the presence of heavy flavors are not without 
intricacies. In particular, care in the choice of a proper factorization scheme is essential for the efficient 
separation of short- and long-distance contributions to the heavy quark cross section. Recent reviews of the 
theoretical status can be found in Refs. |^, ||, [||. The key issue here is whether, for a given renormalization 
and factorization scale fi ~ Q, the heavy quarks of the JV-th flavor are treated as partons in the incoming 
JL , proton, i.e., whether one calculates the QCD beta-function using N active flavors and introduces a parton 

distribution function (PDF) for the N-th flavor. A related, but separate, issue is whether the mass of 
the heavy quark can be neglected in the hard cross section and/or parton distribution function without 
ruining the accuracy of the calculation. Currently, several factorization schemes J|, ||, @, ||, |[ [To[ [n], [l^] 
are available that differ in the treatment of these issues; of those, the variable flavor number schemes reduce 
to the conventional massless MS factorization scheme in the limit Q 2 M , The general proof of the 
factorization in the presence of nonzero heavy quark masses was given in Ref. 

In the near future, we expect the quality of the heavy flavor data to greatly improve. By 2006, the upgraded 
ep collider HERA will accumulate an integrated luminosity of 1 fb~ , i.e., more than eight times the final 
integrated luminosity from its previous runs. Studies of heavy quarks in ep collisions are also envisioned at 
the proposed Electron Ion Collider [Q and THERA jTl|. By the end of these experiments, a large amount 
of information will be available both for the inclusive observables (like the charm component of the DIS 
structure functions) and differential distributions. 

In our work Jl(|, we address a theoretical formalism for the description of heavy-quark differential distri- 
butions both near the heavy-quark mass threshold (Q 2 w M 2 ) and far beyond it (Q 2 3> M 2 ). Our attention 
is primarily focused on the angular or transverse momentum distributions of the final-state heavy mesons. 
The behavior of such distributions is quite complex, since three distinct momentum scales (the heavy quark 
mass M, the photon virtuality Q, and the transverse momentum of the heavy quark) are involved. As a 
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Figure 1: The neutral current semi-inclusive production e + p — ► e + H + X of heavy mesons in the 7*f> cm. reference 
frame. g M ,p M , and are momenta of the virtual photon, proton, and the heavy meson, respectively. The polar 
angle 8h of the heavy meson defines the scale q% discussed in the text. The resummation effects considered here are 
important in the limit 9h — * 0, i.e., when the final-state heavy meson closely follows the direction of its escape in the 
lowest-order flavor-excitation process 7* + q — > q. 



result, large logarithms of ratios of two very different scales tend to appear in the fixed-order perturbative 
cross sections in some regions of the phase space. These logarithms break the convergence of the series 
in the strong coupling as- To obtain stable theoretical predictions, they have to be summed through all 
orders of as ■ The "resummation" is realized by the reorganization of factorized cross sections and solution 
of equations for the gauge and renormalization group invariance. 

The inclusive and semi-inclusive heavy-flavor cross sections may contain large logarithms of two kinds. 
The logarithms ln(Q 2 /M 2 ) appear already in the inclusive observables and are resummed in the heavy- 
quark PDFs. The resummation of this type constitutes the essence of the massive variable flavor number 
factorization schemes, as, for instance, the ACOT scheme Logarithms of the form (a 7 g/q'^)ln m (q 2 ^/Q 2 ), 
< m < 2n— 1 appear at Q 2 ^> M 2 , q\ when calculating angular differential distributions. Here qr = pr/z 
is the transverse momentum pr of the heavy meson H rescaled by the final-state fragmentation variable 
z = (p ■ ph)/(p ■ q)', P^,q^, and p H are the momenta of the initial-state proton, virtual photon, and heavy 
meson, respectively. The limit q\ <§C Q 2 is equivalent to the limit Oh —* 0, where Oh is the polar angle of H 
in the j*p cm. reference frame (cf. Fig. |l|). If M is neglected (M 2 -C Q 2 ), qr is simply related to Oh- 



x J 1 + cos Oh \x 

The logarithms \n(q^/Q 2 ) originate from the incomplete cancellation of soft singularities and factorization 
of quasi-collinear singularities in the fixed-order cross sections. 1 The formalism for their summation was 
proposed by Collins, Soper, and Sterman (CSS) to describe angular correlations in e + e~ hadroproduction 
ps| [l9| and transverse momentum distributions in the Drell-Yan process ^o|. Recently this formalism was 
extended to describe semi-inclusive DIS (SIDIS) production of practically massless hadrons 22, 2j|. The 
resummation in massless SIDIS serves as a high-energy limit for our calculation. 



A discussion of integrahle singularities in the fixed-order differential distributions for the semi-inclusive DIS charm production 
can be found in Ref. 1171. 
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Figure 2: Balance of various terms in the ACOT scheme and resummed cross section. Graphs show qr on the a;-axis 
and dajdq\ on the y-axis. 



The key result of the CSS formalism is that all logarithms due to the soft radiation in the high-energy limit 
can be summed into a Sudakov exponent. This result can be summarized in the following master equation: 

WS? T = 7 / (0 e * 5 (C< " f) (C "" 8 d) e " + °'° - ° Asr ' (1) 

Here b is the impact parameter conjugate to qr, f and d are parton distributions and fragmentation functions, 
respectively. C m , C out contain perturbative corrections to contributions from the incoming and outgoing 
hadronic jets, respectively. The factor e~ s is the Sudakov exponential, which includes an all-order sum of 
perturbative logarithms ln m (q T /Q 2 ) at b < 1 GeV -1 and nonperturbative contributions at b > 1 GeV" 1 . 
Finally, o~fo is the fixed-order expression for the considered cross section, while oasy (asymptotic piece) 
is the perturbative expansion of the 6-space integral up to the same order of as as in apo- At small qr, 
where terms ln m (q T /Q 2 ) are large, ofo cancels well with uasy, so that the cross section (Q) is approximated 
well by the 6-space integral. At qx > Q, where the logarithms are no longer dominant, the 6-space integral 
cancels with ctasy, so that the cross section (|l|) is equal to opo up to higher-order corrections. 

To extend the CSS resummation from the massless to the massive case, the logarithms of both types 
In (Q 2 /M 2 ) and ln(g^/Q 2 ) should be resummed simultaneously. For the resummation of the logarithms 
ln(<3 2 /Af 2 ), we use the simplified ACOT factorization scheme fl3"| . With regards to the logarithms 
In (qp/Q 2 ) , we observe that the proof of the resummation formula in the original papers ||l8| , |l9| |2(J does not 
rely on the dimensional regularization for the soft and collinear singularities appearing in individual Feynman 
diagrams. Although the dimensional regularization does provide a convenient way to regulate ultraviolet 
singularities in as and PDFs (and hence, defines the MS regularization and factorization schemes), it does 
not have to be invoked to regulate the singularities in the infrared region: other methods (e.g., regularization 
with nonzero quark and gluon masses) will do just fine. As a result, the 6-space integral is well-defined both 
for nonzero quark masses and in the limit M/Q — » 0; hence it cannot contain negative powers or logarithms 
of M with the exception of the collinear logarithms resummed in the PDFs or fragmentation functions. Also 
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Figure 3: The angular distributions of B mesons in the 7*p cm. frame close to the mass threshold (Q = 5 GeV, left 
figure) and at Q 2 S> M 2 (Q = 50 GeV, right figure). In both cases, the cross sections are calculated for x — 0.05 and 
M = 4.5 GeV using the CTEQ5HQ PDFs f|| and Peterson fragmentation functions |§f] with e = 0.0033. The plots 
show the fixed-order cross section (long-dashed line), the 6-space integral (dot-dashed line), the asymptotic piece 
(dotted line), and the full resummed cross section (solid line). 



the quark masses regulate collinear singularities, so that at least some of the logarithms are not dominant 
in the threshold region. The remaining issue is about the accuracy of the approximation. The evolution 
equations used in the derivation of the 6-space integral in Refs. Jl^, did not keep track of the power- 
suppressed terms 0((M/Q) n ), n > 0. As a result, the 6-space integral does not necessarily approximate the 
small-gT behavior of the cross section if M/Q is not negligible. However, the full Eq. (0) does a better job 
and approximates the cross section of the order 0{a k s ) up to corrections of the order 0(a^ +1 , M/Q). This 
happens because terms up to the order 0(a|, M/Q) are canceled between the 6-space integral and ctasy', 
the terms ln m (q?/Q 2 ) that were not resummed (if any) cancel between gfo and ctasy ■ 

Fig. H qualitatively illustrates the balance of various terms in the resummation formula in various regions 
of phase space for the lowest-order contributions. First, consider a fixed-order calculation in the simplified 
ACOT scheme near the threshold. In this region, the cross section da/dq^ is well approximated by the 
0(as) photon-gluon diagram (Fig.||a). To this diagram, we add the lowest-order -f*Q term, which resums 
powers of ln(Q 2 /A/ 2 ) and contributes at qr = (Fig.||b). We also subtract the overlap between the two 
diagrams (Fig. ||c), which is approximately equal to, but not the same as, the j*q contribution. The resulting 
distribution (Fig.||d) is close to the fixed-order result, but has discontinuities in the small qr region. These 
discontinuities are amplified when Q increases. 

In the gT-resummed cross section, the fusion diagram still dominates near the threshold, but now the 
resummed cross section and subtracted asymptotic piece are smooth functions, which cancel well at all 
values of qr (Fig.||e-h). Hence, the distribution is physical in the whole range of qr- Finally, at Q 2 ^> M 2 
the small-^T region the j*g fusion contribution is dominated by the term (Fig.^i). Such singular terms 
are summed through all orders in the 6-space integral corresponding to Fig.|^j and are canceled in the fusion 
contribution by subtracting the asymptotic piece (Fig.||k). 

To further simplify the calculation, we use the results of Refs. Jl^, |||, which to neglect masses of heavy 
quarks entering the hard scattering subprocess directly from the proton. This approximation, which differs 
from the complete mass-dependent cross section by higher-order terms, significantly simplifies the analysis. 
We therefore drop heavy-quark masses in hard scattering subdiagrams with incoming heavy quark lines. 
This approximation leads to massless expressions for the perturbative Sudakov factor and C-functions for 
quark-initiated subprocesses, and mass-dependent C-functions for gluon-initiated subprocesses. 
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Figure |3| demonstrates how various terms in Eq. (|l|) are balanced in an actual numerical calculation. This 
Figure shows the bottom quark production cross section vs. the polar angle Oh of the bottom quark in the 
-f*p cm. frame for the lowest-order processes 7*6 —> b and 7*5 — > bb. At this order, the perturbative Sudakov 
factor is identically zero. Due to the heavy mass of the bottom quark, the 6-space integral is practically 
insensitive to contributions from the nonperturbative QCD region b > 1 GeV -1 , so that it can be calculated 
without introducing a phenomenological parameterization for the Sudakov factor in the nonperturbative 
region. The same calculation can be done for the charm production, but in that case the nonperturbative 
Sudakov factor cannot be neglected. The left figure shows that near the heavy-quark mass threshold (Q = 5 
GeV) the 6-space integral cancels well with its perturbative expansion <jasy , so that the the full cross section 
is practically indistinguishable from the fixed-order term. On the other hand, at Q = 50 GeV (right figure) 
the full cross section is substantially larger than the fixed-order term at Oh < 10° and is dominated by the 
6-space integral. In this region, ofo is canceled well by oasy ■ At Oh > 10° the cross section agrees well 
with the massless result. More details can be found in our upcoming publication. 

To summarize, we presented a method for the calculation of angular distributions in the heavy quark 
production both near the heavy-quark threshold and at large momentum transfers. Our method is realized 
in a massive variable flavor number factorization scheme and resums transverse momentum logarithms with 
the help of the 6-space resummation formalism ||l8| , fts| , po|| . Our approach provides an adequate tool for 
detailed studies of differential distributions in a wide range of Q that will soon be available from HERA. 

This work is supported by the U.S. Department of Energy, the National Science Foundation, and the 
Lightner-Sams Foundation. 
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